Comparative study of the electronic structure, phonon spectra and electron-phonon 

interaction of ZrB 2 and TiB 2 
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The electronic structure, optical and x-ray absorption spectra, angle dependence of the cyclotron 
masses and extremal cross sections of the Fermi surface, phonon spectra, electron-phonon Eliashberg 
and transport spectral functions, temperature dependence of electrical resistivity of the ME>2 (M=Ti 
and Zr) diborides were investigated from first principles using the full potential linear muffin-tin 
orbital method. The calculations of the dynamic matrix were carried out within the framework of 
the linear response theory. A good agreement with experimental data of optical and x-ray absorption 
spectra, phonon spectra, electron-phonon spectral functions, electrical resistivity, cyclotron masses 
and extremal cross sections of the Fermi surface was achieved. 
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I. INTRODUCTION 



Ceramics based on transition metal borides, nitrides, 
and carbides have extremely high melting points (>2500 
°C) and are referred to as ultra-high temperature 
ceramics^ 2 - Among them, diborides such as ZrB2 and 
H1B2 have a unique combination of mechanical and phys- 
ical properties: high melting points (>3000 °C); high 
thermal and electrical conductivity; chemical inertness 
against molten metals; great thermal shock resistance^— 
Thus, although carbides typically have the highest melt- 
ing points (>3500 °C), the diborides ZrB2 and HfB2 are 
more attractive candidates for high-temperature ther- 
momechanical structural applications at temperatures 
>3000 °C»ii2 Potential applications include thermal pro- 
tective structures for leading-edge parts on hypersonic 
re-entry space vehicles,—^ propulsion systems r& fur- 
nace elements,— refractory crucibles, 5 and plasma-arc 
electrodes^ In particular, ZrP>2 has the lowest theo- 
retical density among the ultra-high temperature ceram- 
ics, which makes it an attractive material for aerospace 
applications^ 2 -^ Titanium diboride is also potentially 
useful because it has many interesting physical proper- 
ties, such as low density and unusual strength. 7 T1B2 is 
widely accepted for applications including microelectron- 
ics, diffusion barriers, wear- and erosion-resistant coat- 
ings for cutting tools and other mechanical components. 
In these applications, the material's high hardness, high 
melting point, good electrical conductivity, and acid and 
radiation stability is exploited. 8 

The discovery of superconductivity in MgB2 at 39 K 
by Akimitsu^ has lead to booming activity in the physics 
community and activated a search for superconductivity 
in other diborides. Natural candidates for this search are 
AB2-type light metal diborides (A = Li, Be, Al). How- 
ever, up to now superconductivity has not been reported 
in the majority of these compounds. 10 Only very recently 
has superconductivity below 1 K (T c = 0.72 K) been re- 
ported in BeB 2 .75.— According to Ref. [H no supercon- 



ducting transition down to 0.42 K has been observed in 
powders of diborides of transition metals (A = Ti, Zr, Hf, 
V, Ta, Cr, Mo, U). Only NbB 2 is expected to supercon- 
duct with a rather low transition temperature (< 1 K), 
and contradictory reports about superconductivity up to 
T c =9.5 K in TaB 2 can be found in Ref. OJ. Finally, the 
reported T c —7 K in ZrB^ encourages further studies of 
these diborides. 

Presently, a number of experimental studies ex- 
ist dealing with the physical properties of ZrB2 and 
TiB2 such as electric transport properties j 7 -^ - — the 
dc Haas-van Alphen (dHvA) measurements of the 
Fermi surfacejiSr— optical^ 2 - - — and electron-energy- 
loss spectra , 25 ' 26 x-ray absorption and photoemis- 
sion spectra ) 27 ' 28 magnetic susceptibility 2 ?^ and NMR 
measurements^ the phonon density of states^ and 
electron-phonon interaction^ 2 -^ 3 - First-principles calcula- 
tions of the electronic structure of diborides have been 
widely presented ^—^S^— 

Despite a lot of publications, there are still many open 
questions related to the electronic structure and phys- 
ical properties of transition metal diborides. Ihara et 
al£L calculated the band structure and the density of 
states (DOS) of ZrB2 by using an augmented plane wave 
method. They pointed out that the band structure of 
ZrB2 is determined by the sp 2 hybridization, p z state of 
B and the 4d and 5s states of Zr. Similarly Johnson 3 ^ cal- 
culated the band structure of ZrB2 using the Korringa- 
Kohn-Rostoker method in the spherical muffin-tin ap- 
proximation. However, they concluded that the B 2s 
states are localized and do not hybridize with B 2p. 
Pablo et alS^ compared the electronic structure of iso- 
structural alkaline-earth diborides using a full-potential 
linearized augmented plane wave (FLAPW) method and 
found that Zr-B bonds have covalent character, yet still 
remain highly ionic. Fermi surfaces and DOS values at 
the Fermi level reported by Shein et a?4S and those by 
Rosner et al£^ are quite different. Vajeeston et al£L 
also investigated the electronic structure of AlB2-type di- 
borides using the tight-bonding linear muffin-tin orbital 
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(TB-LMTO) method, they claimed that metal-metal and 
metal-boron interactions are less significant than the p— p 
covalent interaction of boron atoms. Burdett et al.^ 
on the basis of orbital overlap, indicated the importance 
of the interaction of orbitals of the metal with those of 
a graphite-like net of boron atoms as well as the in- 
teraction with those of other metals in influencing the 
properties of these species. The bonding nature, elas- 
tic property and hardness were investigated by Zhang 
et al^ for ZrB2 using the plane-wave pseudopotential 
method. The stiffness and the thermal expansion coeffi- 
cient of ZrB2 were calculated using the density functional 
theory formalism by Milman et al. in Ref. |52| . Kaur 
et al^ studied the cohesive and thermal properties of 
these compounds using the rigid ion model. The elastic 
properties, electronic structure, electronic charge distri- 
bution, and equation of states of titanium diboride were 
studied by Milman and Warren^ Perottoni et al.^ and 
Camp et al. 56 using the first-principles methods. Peng et 
al£^ investigated the thermodynamic properties of TiB2 
using a plane-wave pseudopotential method. Munro^ ex- 
amined the physical, mechanical, and thermal properties 
of polycrystalline TiB2 and showed that these proper- 
ties are significantly related to the density and grain size 
of the used specimens. Deligoz et a/j 48 i 50 investigated 
the structural and lattice dynamical properties of TiB2 
and ZrB 2 together with VB 2 , ScB 2 , NbB 2 , and MoB 2 . 
They specifically presented following properties: lattice 
parameters; bond lengths; phonon dispersion curves and 
corresponding density of states; some thermodynamic 
quantities such as internal energy, entropy, heat capacity, 
and their temperature-dependent behaviors. Systematic 
trends in lattice constants and heats of formation for 
these compounds were studied by Oguchi.— Vajeeston 
et alM- investigated the electronic structure and ground 
state properties of these diborides using TB-LMTO. X- 
ray absorption and photoemission spectra of ZrBo_ and 
TiB2 were measured experimentally in Refs. I27II28II59I . 
The optical spectra of ZrB2 was investigated experimen- 
tally by several authors ^2r— however, there is neither 
theoretical investigation of the x-ray absorption spectra 
or the optical properties of the transition metal diborides. 

The band structure and Fermi surface parameters 
were studied by Shein and ivanovskii 4 ^ using the self- 
consistent full potential linearized muffin-tin orbital (FP- 
LMTO) method for ZrB 2 and NbB 2 . Rosner et alM& 
provided a comparison of full potential band calculations 
of the Fermi surfaces areas and masses of MgB2 and ZrB2 
with dHvA data for several symmetry points in the Bril- 
louin zone (BZ). They found, with one possible excep- 
tion, that LDA provides a good description for ZrB2. For 
MgB2 some disagreement in FS areas can be accounted 
for by a shift of tt (B p z ) bands with respect to a (B 
sPxPy) bands by 240 meV and by a readjustment of the 
"Fermi energies" of each of these bands by ±120 meV. 
Heid et al^ measured the phonon density of states of 
MB 2 with M=Ti, V, Ta, Nb, and Y using inelastic neu- 
tron scattering. Experimental data were compared with 



ab initio density-functional calculations using the mixed 
basis pseudopotential method. The results do not exhibit 
indications of strong electron-phonon interaction in the 
diborides considered. Singh 42 studied electron-phonon 
interaction in ZrB2 and TaB2 using a FP-LMTO method. 
The results for phonon density of states and Eliashberg 
function show electron-phonon coupling in ZrB2 to be 
much weaker than in TaB2. The average electron-phonon 
coupling constant A is found to be 0.15 for ZrB2 and 0.73 
for TaB2. Solutions of the isotropic Eliashberg gap equa- 
tion indicate no superconductivity for ZrB2- 

The aim of this work is a complex comparative in- 
vestigation of the electronic structure, optical and x-ray 
absorption spectra, angle dependence of the cyclotron 
masses and extremal cross sections of the Fermi surface, 
phonon spectra, electron-phonon interaction and electri- 
cal resistivity of the diborides TiB2 and ZrB2. The paper 
is organized as follows. Section HT1 presents the details of 
the calculations. Section llTTI is devoted to the electronic 
structure as well as optical and x-ray absorption spec- 
tra, angle dependence of the cyclotron masses and ex- 
tremal cross sections of the Fermi surface, phonon spec- 
tra, electron-phonon interaction and electrical resistivity 
using the FP-LMTO band structure method. The results 
are compared with available experimental data. Finally, 
the results are summarized in Sec. IIVI 



II. COMPUTATIONAL DETAILS 

Most known transition-metal (M) diborides MB2 are 
formed by group III- VI transition elements (Sc, Ti, Zr, 
Hf, V, Nb, and others) and have a layered hexagonal C32 
structure of the AlB2-type with the space group symme- 
try PQ/mmm (number 191). It is simply a hexagonal 
lattice in which closely-packed transition metal layers are 
present alternative with graphite-like B layers (Fig. [1}. 
These diborides cannot be exactly layered compounds 
because the inter-layer interaction is strong even though 
the M layers alternate with the B layers in their crys- 
tal structure. The boron atoms lie on the corners of 
hexagons with the three nearest neighbor boron atoms 
in each plane. The M atoms lie directly in the centers of 
each boron hexagon, but midway between adjacent boron 
layers. Each transition metal atom has twelve nearest 
neighbor B atoms and eight nearest neighbor transition 
metal atoms (six are on the metal plane and two out of 
the metal plane). There is one formula unit per primi- 
tive cell and the crystal has simple hexagonal symmetry 
(D6h). By choosing appropriate primitive lattice vectors, 
the atoms are positioned at M (0,0,0), B (|, i, i), and B 
(§' I> 5) ln * ne un it cell. The distance between M-M is 
equal to c. This structure is quite close packed, and can 
be coped with efficiently and accurately by the atomic 
sphere approximation method. However, for precise cal- 
culation of the phonon spectra and electron-phonon in- 
teraction, a full potential approximation should be used. 

For a crystal where both the fourfold axis and the mag- 
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FIG. 1: (Color online) Schematic representation of the ZrEb 
structure (upper panel) and top view of Zr (large blue circles) 
and B (small yellow circles) planes in ZrB2 (lower panel). 



netization M are perpendicular to the sample surface, 
and the z-axis is chosen to be parallel to them, the di- 
electric tensor is composed of the diagonal e xx and e zz , 
and the off-diagonal e xy components in the formSS 




(1) 



The various elements i a p are composed of real and 
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* e i% where 



imaginary parts as follows: e a p 
a,j3 = x, y, z; e xx = (n + ik) 2 ; n and k are the refractive 
index and extinction coefficient, respectively. The optical 



conductivity tensor a a p = a^l + icr^a ^ s related to the 



J2) 



dielectric tensor e a p through the equation 
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£ a p{u) = 5 a p H (r a p(w). 



(2) 



The optical conductivity of ZrB2 and TiB2 has been 
computed from the energy band structure by means of 
the Kubo-Greenwood^i linear-response expression: 62 



cr a p(u) 
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where /(e„k) is the Fermi function; huj nn >(k) = e„k — e n 'k 
is the energy difference of Kohn-Sham energies; 7 is the 
lifetime parameter, describing the finite lifetime of the 
excited Bloch electron states; n^ n , are the dipole optical 
transition matrix elements^ A detailed description of 
the optical matrix elements is given in Refs. HH. The 
absorptive part of the optical conductivity was calculated 
in a wide energy range. The Kramers-Kronig transfor- 
mation was then used to calculate the dispersive parts 
of the optical conductivity from the absorptive part. We 
used the value 7=0.6 eV for the interband relaxation pa- 
rameter. 

Within the one-particle approximation, the absorption 
coefficient for incident x-ray of polarization A and 

photon energy hcu can be determined as the probabil- 
ity of electronic transitions from initial core states with 
the total angular momentum j to final unoccupied Bloch 
states 



mIM = DZ)K*»k|nA|*i mj )| 2 «(^-^ -M 

(4) 



rrij nk 

x9(E nk -E F ) 



where $ jmj and Ej m . are the wave function and the en- 
ergy of a core state with the projection of the total angu- 
lar momentum mf, f „k and E n ^ are the wave function 
and the energy of a valence state in the n-th band with 
the wave vector k; ep is the Fermi energy. 

Tl\ is the electron-photon interaction operator in the 
dipole approximation 



IT 



(5) 



a/3 



'aP 



where a are the Dirac matrices, is the A polariza- 
tion unit vector of the photon vector potential, with 
a± = 1/V2(l,±i,0),a|| = (0,0,1). Here, + and - de- 
notes, respectively, the left and right circular photon po- 
larizations with respect to the magnetization direction 
in the solid. Then, x-ray magnetic circular and linear 
dichroism are given by /U+ — //_ and fiu — (fi + + /i_)/2, 
respectively. More detailed expressions of the matrix 
elements for the spin-polarized fully relativistic LMTO 
method may be found in Refs. Ir33l|64. 
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The Eliashberg function (the spectral function of the 
electron-phonon interaction) expressed in terms of the 
phonon linewidths 7 qz , has the form 6 -^ 

a2F{cj) = 9^wt-) £ ^ S{UJ - ^ (6) 

2wJM{e F ) ^ w q „ 

The line-widths characterize the partial contribution 
of each phonon: 



7 qu = 27ro; qi/ l5k+qj',kj| 2<5 ( e jk ~ £F)^(e k +qj' - e F ). 

(7) 

The clcctron-phonon interaction constant is defined as: 



/>OC 7 

A e - ph = 2 / — a 2 F( W ), (8) 

It can also be expressed in terms of the phonons line- 
widths: 



K- P h - \] M , q "-. o ; (9) 

were N(e F ) is the electron density of states per atom 
and per spin on the Fermi level (e F ) and fl'k+qi'kj is the 
electron-phonon interaction matrix element. The double 
summation over Fermi surface in Eq.([7]) was carried out 
on dense mesh (793 point in the irreducible part of the 
BZ) 

Calculations of the electronic structure and physical 
properties of the TiB2 and ZrB2 diborides were per- 
formed using a scalar relativistic FP-LMTO method 66 
with the experimentally observed lattice constants: 
a=3.167A andc=3.529A forZrB 2 ; 67 a=3.03A c=3.229 
A, for TiB 2 . 68 For the LMTO calculations we used 
the Perdew-Wang 69 parameterization of the exchange- 
correlation potential in general gradient approximation. 
BZ integrations were performed using the improved tetra- 
hedron method. 70 Phonon spectra and electron-phonon 
matrix elements were calculated for 50 points in the irre- 
ducible part of the BZ using the linear response scheme 
developed by Savrasovi££ The 3s and 3p semi-core states 
of TiB2 were treated as valence states in separate energy 
windows (for ZrB2 : 4s and Ap). Variations in charge 
density and potential were expanded in spherical har- 
monics inside the MT sphere as well as 2894 plane waves 
in the interstitial area with 88.57 Ry cut-off energy for 
ZrB 2 and 97.94 Ry cut-off energy for TiB 2 . As for the 
area inside the MT spheres, we used 3k— spd LMTO basis 
set energy (-0.1, -1, -2.5 Ry) with one-center expansions 
inside the MT-spheres performed up to l max = 6. 
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FIG. 2: (Color online) Energy band structure and total DOS 
[in states/(cell eV)] of ZrB2 (full blue lines) and TiB2 (dashed 
red lines). 
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FIG. 3: (Color online) Partial DOSs [in states/(atom eV)] of 
ZrB 2 . 



III. RESULTS AND DISCUSSION 

A. Energy band structure 

Figure [5] presents the energy band structure and to- 
tal density of states (DOS) of ZrB2 (full lines) and 
TiB 2 (dashed lines). The partial DOSs ZrB 2 are 
shown in Fig. [3] Our results for the electronic struc- 
ture of ZrB2 and TiB2 are in agreement with earlier 
calculations i 29 ' 30 ' 42 1 48 ~ — A common feature for all transi- 
tion metal diborides is the deep DOS minimum (pseudo- 
gap) at the Fermi energy separating the valence band 
and the conduction band. According to Pasturel et al.^ 
a pseudo-gap arises because of a strong chemical inter- 
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FIG. 4: (Color online) Comparison of the total DOS 
(full line) with photoemission 2 ^. (open circles) and x-ray 
photoemission 27 (full circles) of ZrB2. 
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FIG. 5: (Color online) The theoretically calculated and exper- 
imentally measured 2 * x-ray absorption spectra at B K edge 
of ZrB 2 . 



action. The M-B covalent bonding is believed to be re- 
sponsible for this effect. Fig. [5] includes a comparison 
of the total DOS for ZrB2 and TiB2- In both systems, 
we observe a deep minimum in the DOS at the Fermi 
energy, although the gap appears slightly broader in the 
case of ZrB2- The Zr Ad states in ZrB2 are the dom- 
inant features in the interval from —12.5 eV to 9 eV. 
These tightly bound states show overlap with B 2p and, 
to a lesser extent, with B 2s states both above and be- 
low ep , implying considerable covalency. Higher-energy 
states between 9 eV and 20 eV above e p appear to arise 
from Zr bp and 6s states hybridized with B 2p states. 
The crystal field at the Zr site (D6h point symmetry) 
causes the splitting of Zr d orbitals into a singlet a\ g 
(d 3z 2_i) and two doublets e\ g (d yz and d xz ) and ei g (d xy 
and 4r 2 -j/ 2 )- The crystal field at the B site (D3/i point 
symmetry) causes the splitting of B p orbitals into a sin- 
glet a4 (p z ) and a doublet e2 (p x and p y ). B s states 
occupy a bottom of valence band between —13.1 eV and 
—3.0 eV and hybridize strongly with B p x and p y and Zr 
d yz and d xz states located at —12.5 eV to —0.5 eV. B 
p x andpy states are located between —12.5 eV and —0.5 
eV. B p z states occupied a smaller energy interval from 
—7.5 eV to —0.5 eV with a very strong and narrow peak 
structure at around —4 eV. 



B. X-ray absorption and photoemission spectra 

Experimentally the electronic structure of ZrB2 and 
TiB2 has been investigated by means of photoemission 
spectroscopy ) 27 ' 28 point contact spectroscopy^ x-ray ab- 
sorption spectroscopy ) 28 ' 59 and optical spectroscopy^ - — 

Figure U shows the experimentally measured photoe- 
mission (PES) 28 and x-ray photoemission (XPS) 27 spec- 
tra of ZrB2 compared with the calculated energy dis- 
tribution of total DOS. The calculated DOS has been 
broadened to account for life-time effects and for the 



experimental resolution. The characteristic features of 
the XPS are divided into three parts ranging from the 
Fermi energy to —4.5 eV (peak a), —4.5 eV to —8.5 eV 
(peak b), and —8.5 eV to —13 eV (peak c). The low 
energy peak c arises mostly from the B 2s states and 
partly from the low energy peak of Zr Ad xz>yz states (see 
Fig. [3]). The major peak a close to the Fermi energy 
is derived by Zr Ad states. B p states as well as the 
Zr Ad XZtVZ states contribute to the broad peak c located 
from —4.5 eV to —8.5 eV. Agreement between experiment 
and theory in energy position of major fine structures is 
reasonably well. However, peak b is slightly shifted to- 
ward lower energy in the theory, besides, peak a does 
not split into two peaks as observed in the experimental 
XPS spectrum. On the other hand, the experimental PE 
spectrum 2 ^ measured at 325.26 eV (open circles in Fig. 
S]) has a single peak a in close agreement with the theoret- 
ically calculated DOS. Intensity of the low energy part of 
the PE spectrum is significantly increased due to inelasti- 
cally scattered electrons. The corresponding background 
was extracted from the experimental XPS spectrum^ It 
is interesting to note that the position of the peak b in 
DOS is in better agreement with the PE spectrum than 
observed in the case of the XPS spectrum. 

X-ray absorption spectra (XAS) were measured by 
Tsuda et al& at the B K and Zr M 2 , 3 edges of ZrB 2 . Ti 
K XAS were measured by Chu et al£^ The XA spectra 
in metals at the K edge in which the Is core electrons 
are excited to the p states through the dipolar transi- 
tion usually attract only minor interest because p states 
are not the states of influencing magnetic or orbital or- 
der. Recently, however, understanding p states has be- 
come important due to XA spectroscopy using K edges of 
transition metals gaining popularity. The K edge XAS 
is sensitive to electronic structures at neighboring sites 
because of the delocalized nature of the p states. 

Figure [5] presents the theoretically calculated and ex- 
perimentally measured B K XPS spectra (Is —> 2p tran- 
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FIG. 6: (Color online) The theoretically calculated and exper- 
imentally measured^ x-ray absorption spectra at Ti K edge 
of TiB 2 . 
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FIG. 7: (Color online) The theoretically calculated and ex- 
perimentally measured^ x-ray absorption spectra at Zr Mi;i 
edges of ZrB 2 . 



elements^) only 3d 3 / 2 -states occur as final states for M 2 
XAS. For the M 3 XAS, 4d 5/2 -states also contribute. Al- 
though the 3p3/2 — > ids/2 radial matrix elements are only 
slightly smaller than for the 3p 3 / 2 — > 4d 5 / 2 transitions, 
the angular matrix elements strongly suppress the 3p 3 / 2 
—> 4c? 3 / 2 contribution^ Therefore in neglecting the en- 
ergy dependence of the radial matrix elements, the M 2 
and the M 3 spectra can be viewed as a direct mapping 
of the DOS curve for 4d 3 / 2 - and 4ds/ 2 -character, respec- 
tively. 

C. Optical spectra 

The optical spectra of ZrB 2 have been measured by 
several authors^ - — Currently there are no such mea- 
surements for TiB 2 . Fig. [5] shows the theoretically cal- 
culated and experimentally measured optical reflectivity 
spectra R(uj) as well as dielectric constants Si(u>) and 
e 2 (cj) for ZrB 2 . Also presented are theoretically calcu- 
lated R(oj), £i(uj), and £ 2 (w) for TiB 2 . Theory repro- 
duces well peculiarities of ZrB 2 optical spectra. 

We performed decomposition of the calculated e 2 spec- 
trum into the contributions arising from separate inter- 
band transitions and different places of k space. We 
found that the major peak in the e 2 (w) (around 1 eV) 
is mostly determined by the 5 — >■ 6 interband transitions 
along the T— A and A— L symmetry directions (Fig. [9j. 
The shoulder at 2 eV is due to the 4 — !> 5 interband tran- 
sitions around A symmetry point (pink dotted lines in 
Fig. P. 

Although the band structures ZrB 2 and TiB 2 are very 
similar (see Fig. [5]), their optical spectra visibly differ 
from each other. The experimental measurements of the 
optical spectra of TiB 2 are highly desirable. 

D. Fermi surface 



sitions). The agreement between the theory and the ex- 
periment is excellent. The low energy peak around 190 
eV is due to transitions from the Is core level to the 
mostly B p z states (see Fig. [3]) with some amount of 
the p x and p v states. Fine structure at 194.5 eV re- 
flects the corresponding peak between 6 eV and 7 eV 
above the Fermi level (Fig. [3]). Figure [B] presents the- 
oretically calculated and experimentally measured^ Ti 
K XPS spectra in TiB 2 . The agreement between theory 
and experiment is also quite good; except for a second 
major peak around 28 eV that is slightly underestimated 
theoretically. 

Figure [7] presents theoretically calculated and experi- 
mentally measured Zr M 2 3 XPS spectra (3p — > Ad tran- 
sitions). Agreement between theory and experiment is 
again good. Because of the dipole selection rules (apart 
from the 4si/ 2 -states that which have a small contribu- 
tion to the XAS due to relatively small 3p —> 5s matrix 



The magnetoresistancc and Hall effect were measured 
in early investigations of the Fermi surface (FS) of ZrB 2 
in 1966 by Piper^ He showed that ZrB 2 is a compen- 
sated semimetal with an effective concentration of 0.04 
electrons/cell with no open trajectories. In 1978, the 
dHvA effect was observed and investigated by Tanaka et 
al*£ and an attempt was made to interpret the dHvA os- 
cillations on the basis of calculations of the band struc- 
ture of CrB^ using the "rigid band" approximation. 
The model obtained for the FS of ZrB 2 was later con- 
firmed by improved calculations performed using the 
FLAPW method. 38 Recently the Fermi surfaces of ScB 2 , 
ZrB 2 and HfB 2 , were studied by Pluzhnikov et al^- us- 
ing the dHvA effect. Their results for ZrB 2 are similar 
to previous measurements by TanakaJ^ 

Theoretical calculations show a ring-like electron FS 
around the K symmetry point (Fig. I10p and of a wrinkled 
dumbell-like hole FS at the A point (Fig. [11]) in ZrB 2 . 
The electron FS and hole FS have threefold and sixfold 
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FIG. 10: (Color online) The calculated electron sheets ol the 
Fermi surface around K symmetry point from the 6th energy 
band of ZrB2. 
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FIG. 8: (Color online) Theoretically calculated (solid blue 
lines) and experimentally measured (open circles)^ optical 
reflectivity spectra (upper panel) and dielectric constants £2 
(middle panel) and E\ (lower panel) of ZrB2. For the £1 func- 
tion the contributions of different interband transitions are 
presented. 
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FIG. 9: (Color online) Theoretically calculated 4 — ► 5 (pink 
dotted lines) and 5 — > 6 (red full lines) interband transitions 
in the optical conductivity of of ZrB2. 





FIG. 11: The calculated hole sheets of the Fermi surface at 
the A symmetry point from the 5th energy band of ZrB2. 



symmetries, respectively. These are broadly consistent 
with the Fermi surfaces used by Tanaka^ to interpret 
their dHvA data. TiB2 has very similar sheets of its 
Fermi surface. 

Figure [12] shows the calculated cross section areas in 
the plane perpendicular z direction and crossed A sym- 
metry point for hole FS (upper panel) and T point for 
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FIG. 12: (Color online) The calculated cross sections in the 
plane perpendicular z direction and crossed A symmetry point 
(upper panel) and V point (lower panel) for ZrB2 (full red 
curves) and T1B2 (dashed blue curves). The labels are pro- 
vided as used in the text. 



electron FS (lower panel) of Z1B2 and T1B2. It can be 
clearly seen that TiB2 has a smaller FS than Z1-B2. 

Figure [13] represents angular variations of the experi- 
mentally measured dHvA frequencies^ for ZrB2 in com- 
parison with the first-principle calculations for field di- 
rection in the (1010), (1120), and (0001) planes. The ob- 
served frequencies of a, /?, 7, and S oscillations belong to 
electron FS around the K point (see Fig. 4 in Ref. uM . 
The e, v, /i, and C, orbits belong to the hole wrinkled 
dumbbell FS. The a frequencies have four branches at 
the (1010) plane and three branches at the (1120) plane. 
The lower 7 frequencies have two branches in both the 
planes. 

The theory reasonably well reproduces the frequencies 
measured experimentally. However, there are still some 
discrepancies. For high frequencies in the < 0001 > di- 
rection, we found the e and \x branches but were unable to 
obtain v branch. We also discover a new branch tr which 
is not detected experimentally. This branch belongs to 
the electron FS around the K point. It has almost a con- 
stant frequency at the (1010) plane and rapidly drops in 
frequency at the (0001) plane. The theoretically calcu- 
lated £ orbits exist in a wider angle interval than observed 
experimentally. 
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FIG. 13: (Color online) The calculated (open red and blue 
circles for the electron and hole surfaces, respectively) and ex- 
perimentally measured^ (black full squares) angular depen- 
dence of the dHvA oscillation frequencies in the compound 
ZrB 2 . 
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FIG. 14: (Color online) The calculated (open red and blue 
circles for the electron and hole surfaces, respectively) and ex- 
perimentally measured^ (black full squares) angular depen- 
dence of the dHvA oscillation frequencies in the compound 
TiB 2 . 
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FIG. 15: (Color online) The calculated angular dependence 
of the cyclotron masses for the electron Fermi surface (open 
red circles) and the hole Fermi surface (blue open triangles) 
and experimentally measured ones (black full squares) in the 
compound ZrE>2. 



FIG. 16: (Color online) The calculated angular dependence 
of the cyclotron masses for the electron Fermi surface (red 
open circles) and the hole Fermi surface (blue open triangles) 
in the compound TiE>2. 



E. Phonon spectra 



Figure [14] represents an angular variation of exper- 
imentally measured dHvA frequencies^ in Til$2 when 
compared with theoretically calculated frequencies. The 
theoretical calculations quite well reproduce the angle de- 
pendence of the extremal cross sections for low frequency 
orbits 7, a and j3. Similar to ZrB2 we detected theo- 
retically a new branch a in TiE$2 which is not observed 
experimentally. This branch belongs to the electron FS 
around the K point. We also find an additional orbit tt 
at the (1120) plane which is absent in ZrB2 and did not 
detected experimentally. We were not able to find theo- 
retically low frequency 5 oscillations appeared in a small 
angle interval near the < 1120 > direction. For high fre- 
quencies we found the e, y, and £ branches similar to the 
corresponding orbits in ZrB2. However, these orbits have 
not been detected in the dHvA experiment. 20 One of the 
possible reasons for that is the relatively large cyclotron 
masses for these orbits. Figures [TBI and [TBI show the cal- 
culated angular dependence of the cyclotron masses for 
ZrB2 and T1B2, respectively. The cyclotron masses for 
the e, y, and £ orbits in TiB2 are much higher than the 
corresponding orbits in ZrB2. The masses for the low- 
frequency oscillations a, /3, 7 and S are less than 0.2mo 
for ZrB2 and slightly larger in T1B2. 



The unit cell of TB2 (T=Zr, Ti) contains three atoms, 
which give in general case a nine phonon branches. Fig- 
ure [17] shows theoretically calculated phonon density of 
state for ZrB 2 and TiB 2 . The DOS for both ZrB 2 and 
TiB 2 can be separated into three distinct regions. Based 
on our analysis of relative directions of eigenvectors for 
each atom in unit cell, we find that the first region (with 
a peak in phonon DOS at 29 meV in ZrB2 and 37,5 meV 
in TiB2) is dominated by the motion of the transition- 
metal atoms Zr and Ti, respectively. This region belongs 
to the acoustic phonon modes. The shift of the first re- 
gion in the phonon DOS towards lower frequencies for 
ZrB2 in comparison to TiB2 is due to the higher mass 
of Zr. The second wide region (60-80 meV) results from 
the coupled motion of Zr(Ti) and the two B atoms in the 
unit cell. The E\ u , A2 g , B\ g phonon modes (see Table 
[T| lie in this area. The phonon DOS in the third region 
extends from 88 meV to 103 meV in Z1-B2 and from 105 
meV to 115 meV in T1B2. This is due to the movement 
of boron atoms and is expected since boron is lighter 
than transition metal atoms. The covalent character of 
the B-B bonding is also crucial for the high frequency of 
phonons. The in-plane E2 9 mode belongs to this region. 
The second and third regions represent optical phonon 
modes in crystals. The most significant feature in the 
phonon DOS is a gap around 40 to 60 meV for both 
ZrB2 and TiB2 . This gap is a consequence of the large 
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FIG. 17: (Color online) Theoretically calculated phonon den- 
sity of states (full blue lines) for ZrE>2 and TiB2 and experi- 
mentally measured one for TiB2^ (open circles). Dashed red 
line presents the calculated phonon DOS of ZrB2 by Deligoz 
et al. [U- 



TABLE I: Theoretically calculated phonon frequencies (in 
meV) in the T symmetry point for ZrB2 and TiB2 and ex- 
perimentally measured ones for TiBa^ as well as calculated 
phonon frequencies in ZrB2 calculated by Deligoz et al. Ref. 



Compound 


reference 


Eiu 


A 2 g 


Big 




ZrB 2 


our results 
Ref. [50] 


58.70 
60.61 


63.26 
63.49 


71.0 
67.76 


99.70 
98.45 


TiB 2 


our results 
Ref. [32] 


63 .0 
65.5 


63.5 
66.4 


69.1 
70.0 


110.0 
112.8 



mass difference, which leads to decoupling of transition 
metal and boron vibrations. 

The TiB2 phonon DOS was measured using inelastic 
neutron scattering experiments in Ref. |32| . Our results 
are in good agreement with the experiment (see Fig. 1171 
lower panel). The small discrepancy in the positions of 
main peaks for TiB2 does not exceed accuracy of calcu- 
lation. 

Currently, there are no data concerning the experimen- 
tally measured phonon DOS in ZrB2. So we compare 
our results with theoretically calculated phonon DOS by 
Deligoz et al^ (Fig. [T71 upper panel). Calculations of 
these authors were based on the density functional for- 
malism and generalized gradient approximation. They 
used the Perdew-Burke-Ernzerhof functional^ for the 
exchange-correlation energy as it is implemented in the 
SIESTA code*&2i This code calculates the total ener- 
gies and atomic Hellmann-Feynman forces using a lin- 
ear combination of atomic orbitals as the basis set. The 
basis set consists of finite range pseudoatomic orbitals 
of the Sankey-Niklewsky type££ generalized to include 
multiplezeta decays. The interactions between electrons 
and core ions are simulated with the separable Troullier- 
Martinsi^ normconserving pseudopotentials. In other 
words, they used the so-called "frozen phonon" tech- 



nique and built an optimized rhombohedral supercell 
with 36 atoms. This method is inconvenient for calcu- 
lating phonon spectra for small q-points as well as for 
compounds with large number of atoms per unit cell. 
There is a very good agreement between our calculations 
and the results of Deligoz et ai^ in a shape and energy 
position of two first low energy peaks in the phonon DOS. 
There is only a low energy shift of the third peak by ^5 
meV in our calculations in comparison with results of 
Deligoz et al^ (see also Table I). 



F. Electron-phonon interaction 

Figure [T5] shows theoretically calculated Eliashbcrg 
functions for ZrB2 and TiB2. We find no significant 
difference in the shape and energy position of major 
peaks between phonon DOS values and electron-phonon 
coupling functions in these compounds. Therefore, we 
can conclude that electron-phonon Eliashberg function is 
mostly defined by the shape of phonon DOSs in ZrB2 and 
TiB2- There are no regions with unusually high electron- 
phonon interaction and phonon dispersion curves do not 
contain any soft modes which might be indicative of the 
possible superconductivity in these borides. By integrat- 
ing the Eliashberg function using equation ([8]), we esti- 
mate the average electron-phonon interaction constant to 
be Ae-p/^0.14 for ZrB 2 . A similar result was obtained 
earlier by SingbJ2 (A=0.15). Drechsler et a?42 estimated 
the value of the dHvA orbit averaged el-ph coupling con- 
stant to be A <0.1. A weak electron-phonon coupling 
strength of A ^0.1 was derived from both the comparison 
of the calculated density of states at the Fermi level and 
specific heat data (Fuchs et aL— ), and by point-contact 
measurements (Apc=0.06^). 

Figure [TH] (upper panel) represents the PC electron- 
phonon interaction function for ZrB2 in comparison with 
the theoretically calculated Eliashber function. Results 
closely agree in the energy positions of major peaks. 
However, the experimental PC function displays a mono- 
tonically decreasing peak amplitude (as we move along 
the energy scale in the high-energy direction) . As a con- 
sequence, the coupling PC constant Apc"=0.06 is less 
than that obtained from the integration of the Eliash- 
berg function (A e _ p h=0.14). The disagreement might be 
explained by the fact that PC and the Eliashberg func- 
tions have a slightly different nature. First, the kine- 
matic restriction of electron scattering processes in a PC 
is taken into account by a factor K = i(l — ^tan^), 
where 9 is the angle between initial and final momenta 
of scattered electrons (for the Eliashberg function, the 
corresponding factor K=l). Therefore in PC spectra the 
large angle (9 —> n) backscattering processes are dom- 
inated. The second reason for suppressing high-energy 
peaks in the PC function is a deviation from the ballistic 
electron flow in point-contact spectroscopy. (PC spectra 
can not be described in the framework of ballistic regime 
for a high-energy phonon area)^ 
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FIG. 18: (Color online) The theoretically calculated Eliash- 
berg function aF(ui) of ZrB2 and Ti2B2 (full blue lines) and 
experimentally measured point contact spectral function^ 
(open circles) for Zrf^. Dashed red line presents Eliashberg 
function of TiB2 calculated by Heid et al. Ref. [3^ ]. 



FIG. 19: (Color online) Theoretically calculated for the < 
0001 > direction (blue curves) and the basal < 1010 > direc- 
tion (red curves) and experimentally measured temperature 
dependence of electrical resistivity of ZtB^t— (upper panel) 
and TiBz^ (lower panel). 



For TiB2 we again obtain a small value of the electron- 
phonon constant A e - p /i=0.15. Due to absence of the ex- 
perimentally measured electron-phonon spectral function 
in TiB2 we compare our calculations with theoretical re- 
sults obtained by Heid 3 — who used the mixed basis pseu- 
dopotential methodi 80 i 81 There is relatively good agree- 
ment between our calculations and Heid's results for the 
energy position and shape of the peaks (Fig. IT51) . The 
first two major low energy peaks of the Eliashberg func- 
tion are slightly shifted towards the smaller energies in 
comparison with the results of Heidi^ 



G. Electrical resistivity 

In the pure metals (excluding low-temperature region), 
the electron-phonon interaction is the dominant factor 
governing electrical conductivity of the substance. Using 
lowest-order variational approximation, the solution for 
the Boltzmann equation gives the following formula for 
the temperature dependence of pi(T): 



Pi(T) 



irn ceU k B T r°° duj e 2 

N(e F )(vj) J u sinh 2 £ 



c4F( W ), (10) 



where, the subscript / specifies the direction of the elec- 
trical current. In our work, we investigate two direc- 
tion: [0001] (c-axis or z direction) and [1010] (a-axis or 
x-direction) . (vj) is the average square of the / compo- 
nent of the Fermi velocity, £ = Ld/2ksT. 

Mathematically, the transport function at r F(uj) dif- 
fers from aF(w) only by an additional factor [1 — 
w/(k)v/(k')/(w|)], which preferentially weights the 
backscattering processes. 

Formula (fT0|) remains valid in the range 0* r /5 < T < 
26 t r^ where: 



n — / 2\i/2 



(11) 



2 f°° 

{u 2 ) tr = — / we&F(w)dw, (12) 

Mr Jo 



2 / al r F{u) 
Jo w 



(13) 



The low-temperature electrical resistivity is the result 
of electron-electron interaction, size effects, scattering on 
impurities, etc., however, for high temperatures it is nec- 
essarily to take into account the effects of anharmonicity 
and the temperature smearing of the Fermi surface. The 
9 tr =604.8 K and 646.19 K for ZrB 2 and TiB 2 , respec- 
tively. 

Figure [T!J] represents the experimental data for mono- 
crystalline ZrB^ as well as our calculations (upper- 
panel). No evidence of anisotropy of the electrical re- 
sistivity was found experimentally. Our theoretical cal- 
culations also show quite small anisotropically behavior 
of the electrical resistivity in ZrB 2 (compare red and blue 
curves in Fig. 1191 upper panel). There is a good agree- 
ment between our calculations and experimentally mea- 
sured results in the region up to 350 K. 

We found that the anisotropy of the electrical resistiv- 
ity in TiB 2 (Fig. 1191 lower panel) is larger than it was in 
ZrB 2 . Our theoretical results slightly exceed experimen- 
tal data^ especially at high temperatures. This is due to 
using in our calculations the lowest-order variational ap- 
proximation in solution of the Boltzmann equation which 
gives upper limit for the electrical resistivit y 65 ' 83 
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IV. SUMMARY 

We have studied the electronic structure and physical 
properties of ZrB2 and TiE$2 using a full potential linear 
muffin-tin orbital method. We investigated the electron 
and phonon subsystems as well as the electron-phonon 
interaction in these compounds. The theory shows good 
agreement with experimentally measured x-ray absorp- 
tion spectra at the B and Ti K and Zr edges. Agree- 
ment between the experiment and the theory in optical 
spectra of ZrB2 is also good. We found that the major 
peak in the £2^) of ZrB2 around 1 eV is mostly deter- 
mined by the 5 — > 6 interband transitions along T— A and 
A— L symmetry directions. 

We investigated the Fermi surface, angle dependence of 
the cyclotron masses, and extremal cross sections of the 
Fermi surface of ZrB2 and TiB2 in details. Theoretical 
calculations show a ring-like electron FS in ZrB2 around 
the K symmetry point and a wrinkled dumbbell-like hole 
FS at the A point. T1B2 has a smaller FS than ZrB2. 
Theory reproduces the experimentally measured dHvA 
frequencies in both the ZrB2 and TiB2 reasonably well. 
We found that masses for low-frequency oscillations a, /3, 
7, and 5 are less than 0.2mo- Masses for high-frequency 
oscillations e, v, /z, and C are large. We discover new 
branches a both in ZrB2 and T1B2 which did not detected 
experimentally. Theoretical calculations closely repro- 
duce the angle dependence of the extremal cross sections 
of high frequency orbits e, /i, and £ in ZrB2. Similar or- 
bits appeared in the theoretical results for TiB 2 , but not 
detected experimentally. The cyclotron masses for these 
orbits in TiB2 are much higher than the corresponding 
orbits in ZrB2 (compare Figs. (fT5]> and (fTB)) ). It could 
be one of the reasons why they have not been observed 



in the dHvA measurements^ 

Calculated phonon spectra and phonon DOSs for both 
ZrB2 and TiB2 are in good agreement with experimen- 
tal results as well as previous calculations. The Elish- 
berg function of electron-phonon interaction in ZrB2 is in 
good agreement with the experimentally measured point 
contact spectral function for both the position and the 
shape of the major peaks. We did not find regions with 
high electron-phonon interaction or phonon dispersion 
curves with soft modes in either ZrB2 or T1B2. This is 
in agreement with the fact that no trace of superconduc- 
tivity was found in these borides. The averaged electron- 
phonon interaction constant was found to be rather small 
A e - P h=0.14 and 0.15 for ZrB2 and TiB2, respectively. We 
calculated the temperature dependence of the electrical 
resistivity in ZrB2 and TiB2 in the lowest-order varia- 
tional approximation of the Boltzmann equation. We 
found rather small anisotropical behavior of the electri- 
cal resistivity in ZrB2 to be in good agreement with ex- 
perimental observation. We found that the anisotropy of 
electrical resistivity in TiB2 is larger than it is in ZrB2. 
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